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“A Global Magnetohydrodynamic Model of Jovian Magnetosphere” 

The goal of this project was to develop a new global magnetohydrodynamic model of the 
interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits 
of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 
11, Voyager 1 and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, 
complicated system. The Jovian aurora also has been monitored for several years. Like 
auroral observations at Earth, these measurements provide us with a global picture of 
magnetospheric dynamics. Despite this wide range of observations, we have limited 
quantitative understanding of the Jovian magnetosphere and how it interacts with the 
solar wind. For the past several years we have been working toward a quantitative 
understanding of the Jovian magnetosphere and its interaction with the solar wind by 
employing global magnetohydrodynamic simulations to model the magnetosphere. Our 
model has been an explicit MHD code (previously used to model the Earth s 
magnetosphere) to study Jupiter’s magnetosphere [Ogino et al. 1998], We continue to 
obtain important insights with this code, but it suffers from some severe limitations. In 
particular with this code we are limited to considering the region outside of 15Rj, with 
cell sizes of about 1.5Rj. The problem arises because of the presence of widely separated 
time scales throughout the magnetosphere. The numerical stability criterion for explicit 
MHD codes is the CFL limit and is given by At/ Ax < 1 where C m(a is the maximum 
group velocity in a given cell, Ax is the grid spacing and At is the time step. If the 
maximum wave velocity is C w and the flow speed is Cf, C ^ = C w + Cy . Near Jupiter the 

Alfven wave speed becomes very large (it approaches the speed of light at one Jovian 
radius). Operating with this time step makes the calculation essentially intractable. 
Therefore under this funding we have been designing a new MHD model that will be able 
to compute solutions in the wide parameter regime of the Jovian magnetosphere. This 
work was carried out in collaboration with Dr. Jon Linker at SAIC. 

Several years ago we developed the Magnetohydrodynamic Algorithm for Planets (MAP) 
to model the interaction between Jovian plasma and the moons of Jupiter [Kivelson et al., 
1996a,b, Linker et al, 1998, Kivelson et al., 1998]. The form of the MHD equations that 
we solve is given by: 
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Here A is the vector potential, B = V x A is the magnetic field, J = V x B is the current 
density, p is the plasma density, v is the velocity, P is the pressure, y is the ratio of 
specific heats, ti is the resistivity, S is the number of ionization events times the ion mass 
per unit volume per unit time, and v s is the local velocity of neutrals prior to ionization. 
Viscous terms (necessary for capturing sharp fronts) also appear in equations (3) and (4), 
but we have omitted them here for simplicity. The equations are solved in spherical 
coordinates on a non-uniform mesh allowing grid points to be distributed more finely in 
regions of interest. Temporal derivatives are advanced with leapfrog time differencing 
combined with a semi-implicit method (described below); centered spatial differencing is 
used for the terms appearing on the right hand side. 

In the semi-implicit approach we introduce a term into the momentum equation (4) that 
effectively modifies the inertia of the short-wavelength modes while accurately treating 
the long wavelengths, enabling the time step to exceed the CFL limit. The code takes 10 
to 50 times less time to run than a comparable explicit method. When the time step is 
reduced below the CFL limit the method becomes identical to an explicit method. 

Under funding from this grant, along with our partners at SAIC, we developed a 
magnetospheric version of this code. The code has successfully been used to model a 
magnetosphere without corotation at Jupiter’s orbit. The next step in this process will be 
to add rotation to the simulation code. Note while the grant at UCLA has ended the 
companion contract at SAIC has a few months to run. This will enable us to get the 
rotating version working. 

We plan to calibrate our new code by using our older code in the regions where the older 
code is valid. To that end we have used the older code to study the effects of the solar 
wind dynamic pressure and interplanetary magnetic field (IMF) on the configuration of 
the Jovian magnetosphere [ Walker et al., 2001]. We found that both the solar wind 
dynamic pressure and the IMF can cause substantial changes in the magnetosphere. On 
the dayside when the pressure increases the bow shock and magnetopause move toward 
Jupiter and the equatorial magnetic field in the middle magnetosphere becomes more 
dipole-like. When the pressure decreases the boundaries move farther from Jupiter and 
the dayside magnetic field becomes stretched out into a more tail-like configuration. For 
northward IMF the boundaries move toward Jupiter but the field becomes more tail-like 
while for southward IMF the boundaries move away and the field becomes more dipole- 
like. These changes are qualitatively consistent with those observed on spacecraft 
passing through the dayside magnetosphere. However, we were not always able to get 
quantitative agreement. In particular the model does not reproduce the extremely tail-like 
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magnetic field observed during the Pioneer 10 and Ulysses inbound passes. Tailward 
flows were found in the nightside equatorial plasma sheet for most IMF orientations. 
Both inertial effects and the IMF influence reconnection in the tail. The only time the 
tailward flow in the magnetotail stopped was during prolonged intervals with southward 
IMF. Then reconnection in the polar cusp caused the flow to move out of the equatorial 
plane. 

We also began using the Jupiter simulations to interpret auroral observations [Waite et 
al., 2001]. On Septerber 22, 1999 the Hubble Space Telescope observed a rapidly 
evolving, very bright and localized emission poleward of the main auroral oval. We used 
the simulation to model the mapping of the auroral flare along magnetic field lines from 
the polar ionosphere into Jupiter’s magnetic equator and found it was connected to 
Jupiter’s outer magnetosphere. 
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